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Abstract 

Using the frequency-domain theory, the high-order above threshold ionization (HATI) of a 
molecule is treated as a two-step transition: an above threshold ionization (ATI) followed by a 
laser-assisted collision (LAC). It is found that the HATI spectrum carries three pieces of infor- 
mation: the fingerprint of the molecular wave function lying in the ATI transition, the geometry 
structure of the molecule lying in the potential scattering between two plane waves, and the inter- 
action between the ionized electron and the laser field lying in a phase factor in the LAC transition, 
where the phase factor can be regarded as due to the action difference of the classical trajectory 
before and after the collision within one laser cycle. Based on this picture, it is demonstrated that 
the skeleton structure of a molecule can be imaged more clearly using a relatively higher frequency 
laser field due to the simplification of interference originated from the phase factor. This work can 
shed light on imaging the structure of a complex molecule by a HATI process. 
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Laser-driven rescattering process of electrons opens a new way to image the structure 
of a molecule. Both high-order harmonic generation (HHG) and high- 

order above threshold ionization (HATI) processes have been utilized to realize this 

task. Odzak and Milosevic [6j have shown clear minima on HHG spectra that depend on 
molecular orientation. Becker et al. [8] have demonstrated that the interference structure of 
a HATI spectrum strongly depends on the symmetry of the ground state molecular orbital. 
A general formula that describes the destructive interference of a multiple-center recollsion 

n ri 

has been found in both HHG p] and HATI [8j spectra. Now people are trying to extract 
more information about a molecule from the HATI spectrum. Therefore, it is necessary to 
understand in more detail how a complex interference pattern forms in the process, and 
what characteristic of the whole system the pattern represents. 

The frequency-domain theory based on a nonperturbative quantum electrodynamics was 

n 



developed by Guo et al. [10( in 1989 dealing with the above threshold ionization (A 



Then, Gao et al. 11[ extended this method to a study of HHG process. Fu et al. 12 



'I)- 



further explained the origin of the plateau structure on HHG spectrum and demonstrated 
the relationship between the time-domain and frequency- domain theories in strong laser 



physics. More recently, Wang et al. 14| applied this theory to investigate HATI process. 
It should be noticed that the characteristic of the frequency-domain theory is that the 
time-evolution dynamic process of a quantum system in a strong laser field is mapped 
into a frequency space; hence all the rescattering processes can be treated as a two-step 
transition. Especially, the HATI of a molecule in a strong laser field can be treated as an 
ATI followed by a laser-assisted collision (LAC). Consequently, we find in this paper that 
the characteristic of a molecule-laser system can be identified from the HATI spectrum: the 
fingerprint of the molecular wave function lies in the ATI transition, the geometry structure 
of the molecule lies in the potential scattering term between two plane waves, and the 
interaction between the ionized electron and the laser field lies in a phase factor in the 
laser-assisted transition. Furthermore, since these three parts are relatively independent 
with each other, we demonstrate that the molecular skeleton structure can be imaged more 
clearly with a higher frequency laser field because of the simplification of interference due 
to the interaction between the ionized electron and the laser field. Recently, we are aware of 
the work by Reiss on the ATI of an atom in a low frequency laser field using the strong-field 
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approximation 15]. 



The frequency- domain theory of the HATI process has been described in [ijj]. Here we 
summarize the basic elements of this theory and the modifications for a molecular case. 
Atomic units are used throughout unless otherwise stated. We consider a quantized single- 
mode laser field of frequency u. The Hamiltonian for a molecule-laser system is 

H = H + U(r) + V, (1) 

where H = + uN a is the energy operator for a free electron-photon system, 

N a = (a) a + aa Jf )/2 is the photon number operator with a (a^) being the annihilation 
(creation) operator of the laser photon mode, U(r) is the molecular binding potential that 
characterizes the geometry structure of the molecule, and V is the electron-photon inter- 
action V = — ^-A(r) ■ (— *v) + 6 ; in which A(r) = g(eae lk ' T + c.c.) is the vector 
potential, and g = {2ujV e )~ l l 2 with V e being the normalization volume of the field, and e the 
polarization vector of the laser field. 

The time- independent feature of the field- quantized Hamiltonian enables us to treat HATI 
as a genuine scattering process in an isolated system that consists of the photons and the 
molecule. Since the energy is conserved throughout the interaction, the formal scattering 
theory [161 ] can thus be applied. Therefore, the transition matrix element can be written 



as 10] 



Tft =(^f\V\ij l ) + ^ f \U- -—V\1h), (2) 

hif — ti + le 

where the initial state = |$j(r), r^) = $j(r) ® \rii) is the eigenstate of the Hamiltonian 
Ho + U(r) with the associated energy Ei = —Eb + (n» + \)uj. Moreover, $i(r) is the ground- 
state wave function of the molecular electron with the binding energy Eb and |n») is the Fock 
state of the laser mode with the photon number rij. Finally, the final state \ipf) = \^ Pfnf ) 
of energy Ef = E Pfn . is the Volkov state of the quantized-field [h], Q. 

The first and second terms in Eq.(j2J) respectively correspond to the processes of direct 
and rescattering ATI. Thus, can be expressed as T/j = Td + T r . On the one hand, as for 



the direct ATI transition, we have 
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T d ={^ f \V\^ i ) 

= V- 1 ftu>(u p -j)*(p i )J j (£f,r,), (3) 



3 



where u p = U p /uo with U p the ponderomotive energy of an electron in the laser field, j = 
rii — rif, $(pi) is the Fourier transform of the initial wave function $j(r), and Jj(Cf,v) = 
J2m=-ooJ-3-2m((f)Jm(v) is tne generalized Bessel function with ( f = 2^Ju p /u p f ■ e, p/ 
the momentum of the electron in the finial state, and 77 = u p /2. In general, a molecular 



ground-state wave function can 
(LCAO) approximation as in 



je constructed using the linear combination of atomic orbital 



[s, Q, Q. Since Eq. (J3]) contains a molecular wave function, 
the direct ATI spectrum can thus image the characteristic of a molecular wave function. 
As an example, we calculate the angle-resolved direct ATI spectra for O2 and N 2 molecule 
with the molecule axis along the laser polarization and the laser frequency u = 0.114 and 
U p = 2uj. The angle Of is between the final momentum of the ionized electron and the 
molecular axis. We find that for N2 the ATI rate decreases with Of . For O2, the ATI rate 
is zero at Of = 0. It then increases rapidly with 8f, gets its maximum at about 20 degrees, 
and then decreases with Of. The dependence of the ATI rate on the angle agrees well with 



the results by Chen et al. [191] . 



On the other hand, the rescattering ATI transition can be written as 

pini 

x($ pini I V I §i,ni)6{E f - E pini ) 



T LAC T ATI 5(E f - E Pini ) . (4) 



-ITT 

all channels 



To obtain Eq. (TjJ, we used the completeness relation of the Volkov states l^pmi) and 
assumed that the effect of the binding potential U can be neglected when the electron is in the 
continuum. The physics underlying Eq. (j4]) is clear. Specifically, Tati — (^pmi I V \ ^i,rii) 
represents the direct ATI amplitude, where the ground state electron absorbs n\ photons 
from the laser field and ionizes; whereas Tlac = {^p f n f I U \ $ Pim ) represents the amplitude 
of LAC in which the ionized electron absorbs n\—Uf photons from the field during its collision 
with the nucleus, resulting in the change of the canonical momentum of the electron from 
Pi to pf. Therefore, from the frequency- domain viewpoint, the rescattering ATI can be 
described simply as an ATI followed by a LAC with all ATI channels summed up coherently. 
Furthermore, Eq. (j4|) indicates that Tati provides a weight amplitude for Tlac transition 
in each ATI channel. Because the initial molecular wave function only appears in the Tati 
term, it can only influence the amplitude of LAC transition, resulting in that the wave 
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function can affects the HATI amplitude of each channel, rather than the structure of the 
HATI spectrum. This finding can be confirmed by the HATI spectra shown in 8]. 

The ATI transition Tati in Eq. (jlj) has the same form as Eq. with pf and rif being 
replaced by pi and n 1; respectively. Furthermore, the transition matrix element of LAC can 
be written as 

= K" 1 ^(Ci-C/)(P/|f/|Pi) 
= K" 1 ^(p/I^Ipi) 



X 



f T 

/ dt exp[— i{sut + (Ci — Cf) smut)], (5) 
Jo 



where J s is the Bessel function with s = ri\ — rif and £i = 2yu p /ujpx ■ e. To ob- 
tain the last equation, we use the integral representation of the Bessel function J n (z) = 
t^t f_ exp \—i{n9 — z sin 8)]d8. The transition term is (pf \ U | pi) = J dr exp[— i(pf — 
Pi) • r]U(r). Equation ([5]) indicates that the LAC can be understood as the potential 
scattering between two plane waves, as shown by (pf | U \ pi), with a phase factor 
AG = Jq dtexp[— i(soot + (d — Q) sin oot)} reflecting the phase difference within one op- 
tical cycle T = 2tc/uj between the two states before and after the collision. The existence of 
the phase factor is due to the fact that the recolliding electron oscillates by the laser field 
during the recollision process. Thus, the interference pattern of the molecular HATI spec- 
trum can be attributed to the combination of the potential scattering and the phase factor 
A0. In particular, the skeleton structure of the molecule lies in the potential scattering 
term, while the interaction between the ionized electron and the laser field is represented 
by the phase factor term in Eq. ([5]). Since the phase factor is independent of the potential 
scattering term and is only determined by the laser condition, the phase factor just provides 
a "background" of interference fringes. Based on this perspective, one may expect that the 
structure of a molecule can be imaged more clearly by the HATI spectrum if the "back- 
ground" interference pattern from the phase factor can be simplified. Fortunately, this idea 
can be realized when the laser frequency increases, as shown in Fig{T] 

We now focus on the imaging of the geometry structure of a molecule, rather than its wave 
function, by HATI spectrum. We consider a HATI process of H^" in an intense laser field. We 
employ a zero-range potential model for U(r) similar to that in [20], where the zero-range 
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potential is Uo(r) = ^f5(r)J^r, with k = ^2\Eb\- The two-center binding potential of a 

diatomic molecule can thus be written as U(r) = Uo(r — z ) + Uq(t + z ), where — z and 

+z are the positions of the two nuclei. The molecular orientation is along the polarization 

of the laser field and the ionization potential is 29.8 eV that corresponds to the ionization 

potential of the ground state of with the internuclear separation Rq being 2 a.u. The 

wave function of is $i(r) = , 1 [0(r — z ) + 4>{v + z )], where 0(r) is the atomic 

v 2 (i+c) 

wave function corresponding to the potential Z7 (r) and C = f 0(r — z o )0(r + z )<ir is the 
atomic orbital overlap integral. The laser intensity is 4.7 x 10 14 W/cm 2 and the frequency is 
chosen to be three different values. Figure Q] presents the angle-resolved HATI spectra with 
uj = 0.057 ((a) and (d)), 0.086 ((b) and (e)), and 0.114 ((c) and (f)) for the H+ molecule ((a)- 
(c)) and the atomic xenon ((d)-(f)). The ponderomotive energy U p is 16.0a;, 4.7a;, and 2.0a; 
for the case of u = 0.057, 0.086, and 0.114, respectively. In Fig{TJ 9f is the angle between 
the final momentum of the ionized electron and the molecular axis. It shows that the cutoff 
of the kinetic energy spectrum at Of = is about 10U P for all the three frequency cases. 
Comparing the molecular case with the atomic case, one can find that (1) there are common 
interference fringes, which we define as "background" interference due to A© in both the 
molecular and atomic cases, and these fringes reduce as the laser frequency increases. (2) 
There exist two destructive curves for the molecular case in Fig. [Tj (a-c), which predict 
the minimum positions caused by the destructive interference of the two-center collision, as 
mentioned in 8|, |9|]. (3) The two destructive curves in the molecular case do not change with 
frequency, leading to a much clearer destructive pattern on the HATI spectrum for the case 
of higher frequency, as shown in Fig. [1] (c). 

We first investigate the source of the "background" interference stripes and the reason 
why it decreases with frequency. Let us analyze the LAC transition for separated ATI 
channels. By using the saddle point approximation, Eq. ([5]) can be expressed as follows: 



T hA c= -tHp/ I C/ I pi) 



irV e y oo 2 (d — (f) sino;£i 

x cos[so;ti + (Ci — Q) sino;ti — 7r/4], (6) 

where t± is the saddle point time satisfying cos(o;ti) = (p"j — p\) /[4^Ju p uj (pi — p/) • e], and 
s = (p 2 —pf)/(2Lv) is the number of photons absorbed during the collision. We can find that 
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the interference comes mainly from the cosine function 

cosF = cos[so;ti + (Ci — C/) sincuti — vr/4]. (7) 

The destructive interference lines should occur when F = (2n + 1)tt/2. Figure [2] presents the 
angle-resolved HATI spectra of the Xe atom when u = 0.057 ((a) and (b)) and uj = 0.114 
((c) and (d)) for channel 1 ((a) and (c)) and channel 6 ((b) and (d)) respectively. One can 
find that the interference fringes increase with channel order and decrease with frequency. 
The dashed lines in each panel of Fig. [2] is obtained by the formula F = (2n + 1)tt/2. 
For comparison, we only present the curves with 6j £ (jr/2,7r) in Fig. [2j It indicates that 
the dashed curves by Eq. (jTJ) agree well with the quantum calculations for all the cases 
shown in Fig. [2j except that the curves for the small Ef part in Fig. [2] (b) are a little lower 
than the quantum results. This difference may be due to the error from the saddle-point 
approximation. 

n 

As mentioned in [14], the phase factor in Eq. (jSJ) can be regarded as due to the 
classical action difference of the classical trajectory before and after the collision, i.e., 
AG = J Q T dt exp[— iAS(t, pi, p/)] with AS(t,pi,pf) = S(t, pi) — S(t,pf), where the clas- 
sical action is S(t,p) = | J dt[p — A c (t)] 2 with the potential of the corresponding classical 
laser field being A c (£) = ^^/Up[ee~ lujt + c.c). Using the saddle-point approximation, we may 
find that there are two moments t\ and 2n/u) — t\ that the rescattering occurs within one 
optical cycle, which corresponds to two classical trajectories; hence the interference fringes 
shown in Fig. [2] is due to the interference of these two classical trajectories, where the min- 
ima are determined as F = (2n + l)vr/2. Although Eq. (j7j) is so complex that there is 
no direct way to show the relationship between the number of the interference fringes and 
the system condition, the results in Fig. [2] indicate that some information can be derived 
from the interference fringes. Figure [2] shows that the smaller the U p /u and the smaller the 
incoming kinetic energy of the rescattering electron, the fewer the interference fringes. 

Comparing the total HATI spectra in Fig. [1] with the channel HATI spectra in Fig. [2J one 
can see that, for the case of low frequency, the total HATI spectrum is more complex than 
that of channel 1; while for the case of high frequency, the total spectrum agrees with that of 
channel 1. Since the total spectrum is the coherent summation of the contributions from all 
ATI channels, the above results indicate that the number of ATI channels that contribute 
effectively to HATI decreases as uj increases. To confirm this finding, we present in Fig. [3] the 
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HATI spectra with Of = for different ATI channels for the case of u = 0.057 (a) and 0.114 
(b) respectively. Figure [3] shows that the strength of HATI spectrum decreases more rapidly 
with channel order for the case of higher frequency, leading to a further simplification of the 
total HATI spectrum. 

We now investigate the interference pattern of the HATI spectrum for different geometry 
structures of the Hj molecular ion. Using the two-centerpotential of the molecule, the 



potential scattering of the plane wave can be expressed as 



R 

(p/ | U | pi) oc cos(— [p/cosfl/ -picosfli]), (8) 

where Of is the angle between the molecular axis and the final momentum p/ and 0\ is 
the angle between the molecular axis and the momentum pi. Therefore, the condition 
PfCosOf — picos#i = (2n + 1)-k/Rq corresponds to the destructive position in the HATI 
spectrum, which contains the characteristic of the geometry structure of the molecule. As 
an example, Fig. @] presents the HATI spectra of Hj~ with Ro = 2 (a), 4 (b), and 5 (c) 
respectively and u = 0.114. For comparison, the corresponding spectra with u = 0.057 are 
also presented (Fig. [4] (d-f)). One can find that the HATI spectrum for the low frequency is 
hard to clearly image the structure of the molecule because the background interference is too 
complex. These results agree with the work by Becker et al. where the angle-resolved 
HATI spectra shown in Figs. 6 and 7 of that paper are so complex that the destructive 
curves are hard to be identified although these curves are predicted (solid lines) by Eq. (jSJ) . 
In contrast, the molecular structure can be clearly imaged by the HATI spectrum as the 
frequency increases, as shown in Fig.[4](a-c). Furthermore, the pattern of the HATI spectrum 
changes dramatically with Rq, and the destructive position can be well predicted by the 
formula pf cos Of + p\ = ( - 2n i ^ ) 1 - >7r (n = 1 for solid black curve and n = 2 for dashed black 
curve) and pf cos Of — Pi = (n—1 for dot red curve). 

In summary, we investigate the HATI of a molecule in a strong laser field by the frequency- 
domain theory. The HATI can be treated as a two-step transition: an ATI followed by a 
LAC. Based on this viewpoint, It is found that the characteristic of a molecule-laser system 
can be identified from the HATI spectrum: the fingerprint of the molecular wave function 
lies in the ATI transition, the geometry structure of the molecule lies in the LAC transition 
with a phase factor which represents the interaction between the ionized electron and the 
laser field. Furthermore, we demonstrate that the HATI spectrum of high frequency laser 
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Figure 1: (Color online) Angle-resolved HATI spectra with frequency uj = 0.057 ((a) and (d)), 
0.086 ((b) and (e)), and 0.114 ((c) and (f)) for H+ molecule ((a)-(c)) and Xe atom ((d)-(f)). The 
results are plotted in log scale. 

can clearly image the molecular structure of large internuclear distance whereas it is hard 
to be identified for the low frequency laser because of its complex nature of the interference 
pattern. 
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Figure 4: (Color online) Angle-resolved HATI spectra of HjJ" with the oj = 0.114 ((a)-(c)) and 0.057 
((d)-(f)) for the internuclear distance Rq = 2 ((a) and (d)), 4 ((b) and (e)), and 5 ((c) and (f)). 
The results are plotted in log scale. 
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